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Abstract 

The equation of state of hadron resonance gas at finite temperature and baryon density is 
calculated taking into account finite-size effects within the excluded volume model. Contributions 
of known hadrons with masses up to 2 GeV are included in the zero-width approximation. Special 
attention is paid to the role of strange hadrons in the system with zero total strangeness. A density- 
dependent mean field is added to guarantee that the nuclear matter has a saturation point and 
a liquid-gas phase transition. The deconfined phase is described by the bag model with lowest order 
perturbative corrections. The phase transition boundaries are found by using the Gibbs conditions 
with the strangeness neutrality constraint. The sensitivity of the phase diagram to the hadronic 
excluded volume and to the parametrization of the mean-field is investigated. The possibility of 
strangeness-antistrangeness separation in the mixed phase is analyzed. It is demonstrated that 
the peaks in the K/ir and A/ir excitation functions observed at low SPS energies can be explained 
by a nonmonotonous behavior of the strangeness fugacity along the chemical freeze-out line. 

PACS numbers: 21.65.Mn, 21.65. Qr, 24.85. +p, 25.75.Nq 
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I. INTRODUCTION 



Relativistic heavy-ion collisions represent a powerful tool for studying properties of 
strongly interacting matter in the laboratory. The main goal is to explore the phase di- 
agram of such matter and, in particular, to investigate properties of a new phase - the 
quark-gluon plasma (QGP). This can be done only indirectly, i.e. by comparing predictions 
of different theoretical models with experimental data. The fluid-dynamical model is one of 
the most popular models for describing the relativistic heavy-ion collisions. In this model 
one needs the equation of state (EOS) and transport coefficients as input information. The 
QCD lattice calculations can give reliable results on thermodynamic properties of strongly 
interacting matter only for small baryon chemical potentials. In this case a crossover type 



of the deconfinement transition is predicted It seems that this prediction is confirmed 
by recent RHIC experiments (see e.g. Q]). The structure of the phase diagram at high 



baryon densities remains rather uncertain [3J. However, exactly the baryon-rich matter 
attracts the main interest due to the possibility of the first order deconfinement phase tran- 
sition [4, 5j. The future facility for antiproton and ion research (FAIR) at GSI (Darmstadt) 
is especially focused on heavy-ion collisions with bombarding energies 20-40 AGeV, where 
maximal baryon densities of about 10 times the normal nuclear density hq = 0.15 fm~ 3 are 

n 

expected [6|. Some interesting results concerning dense baryon matter have been obtained 
already at low SPS energies by the NA49 Collaboration 3, S]- 

To understand properties of the deconfinement transition, one needs to know accurately 
the EOS of the hadronic as well as the quark-gluon phases. As demonstrated in Ref. |9|, 
very often arbitrarily chosen models for these two phases do not lead to any phase tran- 
sition. For instance, taking an ideal resonance gas for the hadronic phase and the bag 
model for the QGP, one comes to the paradoxical conclusion that the hadronic phase is 



thermodynamically stable at high temperatures [10J. It is clear that at high densities the 
repulsive interaction between hadrons becomes important and should be explicitly taken 
into account. One can introduce such interaction via the vector meson exchange as done in 
the relativistic mean-field models of the Walecka type llj. However, the repulsive vector 
field in these models is proportional to the net baryon density and, therefore, can not solve 
the problem in the case of baryon-free matter. Another possibility is to take into account 
the finite size of hadrons within the Van der Waals approach fl2 |. Different versions of 
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the EOS of strongly interac ting; matter with finite-size corrections have been considered in 



Refs. 
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191 ]. but the sensitivity of the phase diagram to the choice of 



hadronic volumes was not fully investigated. In this paper we address this problem on the 
quantitative level. This allows us to derive a consistent EOS with the deconfinement phase 
transition. We also study the role of the strangeness neutrality constraint, which previously 



has been investigated only qualitatively [20l. l2lfl. 



The paper is organized as follows. In Sect. II we describe thermodynamic properties of a 
hadron resonance gas in the excluded volume approximation. In this section we show that the 
strangeness neutrality constraint leads to a non-monotonic behavior of the K/tt and A/tt 
multiplicity ratios along the freeze-out line in the chemical potential-temperature plane. 
The bag model for a quark-gluon phase is formulated in Sect. III. The phase diagram of 
strongly interacting matter is studied in Sect. IV. In particular, we investigate its sensitivity 
to the choice of hadronic volumes. In Sect. V we introduce the mean-field interaction of 
baryons in order to implement the liquid-gas phase transition in nuclear matter at low 
temperatures. In Sect. VI we summarize our results and outline possible improvements of 
the model. 



II. HADRONIC PHASE 

A. Hadron resonance gas within the excluded volume approximation 

Let us consider a purely hadronic system in the total volume V assuming local thermo- 
dynamic equilibrium and neglecting the isospin, Coulomb and surface effects. Since hadrons 
are composite particles, their finite sizes should be taken into account at high enough densi- 
ties or temperatures. The intuitive way to implement the finite-size effects is to use the Van 
der Waals prescription which is known also as the excluded volume approximation. Within 
such an approach the hadronic system is still regarded as an ideal gas, but in the volume 
reduced by the volume occupied by constituents, 

V' = V-^2viNi. (1) 

i 

Here i>, and iVj are the excluded volume and the number of hadrons of type i, the sum 
runs over all hadronic species i . In the following we assume that all hadrons have the same 
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radius r^. Then one can estimate the excluded volume per particle as 1/2 of a spherical 
volume with the radius 2rh [ijl 

167T □ ,„ s 

Vi = v = —r s h . (2) 

This approximation should be good enough at particle densities much smaller than the 
density of close packing The quantity v will be considered below as a model parameter. 



According to the analysis of hadronic yields observed in relativistic heavy-ion collisions [18] , 
the reasonable interval of hadronic radii is ~ (0.3 — 0.6) fm. This corresponds to the 
excluded volume range v ~ (0.5 — 3) fm 3 . 

Making the replacement V — > V in the canonical partition sum, one obtains the following 



equation for the pressure of hadronic system 16j 



P = Y,Pi<Jk,T) , (3) 

i 

where Pj(/2i,T) is the partial pressure of ideal gas of z-th hadrons at temperature T and 
chemical potential /Tj . The relation connecting /Tj with the real chemical potential /ij reads 

jIi = Hi-vP. (4) 

In fact, Eqs. (J3j) — (SJ) give the integral equation for P at a given temperature and set of 
chemical potentials {/ij} . In the limiting case v — > one obtains the Dalton law for a mix- 
ture of ideal gases of different hadrons. 

The condition of chemical equilibrium with respect to strong interactions leads to the 
relations 

Hi = BiH + SiHs, (5) 

where Bi = 0,±1 and Si = 0, ±1,±2... are, respectively, the baryon and strangeness 
quantum numbers of the i-th hadronic species, \x and /j,g denote the baryon and strange 
chemical potentials. Therefore, pressure of a chemically equilibrated system is a function of 
only three independent variables /i, T: 

P = P(jji,fji S ,T). (6) 



This function satisfies the thermodynamic relation 12j 



dP = sdT + ndfi + risdfis , (7) 
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where s is the entropy density, n and ns are, respectively, the baryon and strangeness number 
densities. In this paper we consider only the hadronic matter with zero net strangeness. Such 
matter may be created e.g. in collisions of nonstrange projectile and target nuclei. Since 
strangeness is conserved in strong interactions, the produced system must obey the condition 
of strangeness neutrality 

n s (ii,iis,T) = 0. (8) 

In a baryon-asymmetric matter (fi ^ 0) this condition can be satisfied only at nonzero 
strange chemical potential fis ■ It easy to show that fis should be positive in the baryon-rich 
matter. This follows from the fact that the negative strangeness is carried by hyperons while 
the positive strangeness is associated mostly with kaons. 

From Eqs. (EJ) — (JSJ) , and (J7J), one can obtain explicit relations for thermodynamic quanti- 



ties s, n, ns 



18|: 



s= QT=rJ2si&i, T ), (9) 

i 

dP 

n = y = r) j B i n i (Jii,T), (10) 

i 

ns=^r~ = r , y]3 i ni(Jii,T), (11) 

where Sj = dPi/dT and Hi = dPi/dJli are, respectively, the entropy and particle number 
density of the ideal gas of i-th hadrons. The reduction factor r in r.h.s of Eqs. fl9|)- f|TT|) is 
defined as 

r = (1 +v^n t )-\ (12) 

i 

The number density of the i-th species satisfies the relations 

dP „ 1 
Hi = —— = mi < - . (13) 

Ofli v 

By using Eqs. (E])-©, (j9l)- (IT2"l) and the Gibbs-Duhem equation 

e = -P + Ts + /in + fi s n s , (14) 
one obtains the expression for energy density 

e = r ^2ei(jU,T) , (15) 



5 



where q = — Pi + Tsj + /Xjnj is energy density of the ideal gas of z-th hadrons. Below we take 
into account the contributions of stable hadrons as well as mesonic and baryonic resonances 
in the zero-width approximation. In the grand canonical ensemble the thermodynamic 
functions of the relativistic ideal gas of species i are given by the following explicit rela- 
tions (H = c = 1) 



Pi 
n 



2tt 2 



de- 



exp 



± 1 



■(e 2 -m?) 



\ 



(16) 



/ 



where m 8 is the mass of the i-th hadron and ^ is its spin-isospin degeneracy factor. The 
lower sign in r.h.s. of (TToT) corresponds to mesons [B i — 0) [57| and the upper one, to baryons 
(Sj = 1) or antibaryons = — 1). The integrals in Eq. ffTBT) are calculated numerically as 
explained in Appendix. 

In the present work we include contributions of the lightest hadrons with masses 
mi < 2 GeV. Altogether we take into account 59 mesonic and 41 baryonic species listed 
in Ref. j^. This corresponds to 307 different isospin states of mesons, baryons and an- 
tibaryons. Characteristics of hadrons included in our calculations are given in Table [B For 
hadrons with nonzero isospin we use the isospin averaged mass values. Note, that unless 
otherwise stated, we do not include a very broad scalar meson resonance /o(600) with mass 



m ~ 0.6 GeV and width T > 0.6 GeV 

numbers of pions {dj) and kaons (df) produced in decays of the i-t 
bers are calculated by using the decay branching ratios from Ref. 



58j | . The last two columns of Table [J give the average 

i hadron. These num- 



22j[59|. For particles, 
df = 0. In decays of 



stable with respect to strong interactions, we use the values dj 
heavy hyperons only .fT-mesons, not kaons, are produced. In these cases the last column of 
Table [I] gives the contributions of corresponding antihyperons (see entries with asterisks). 
For example, in the case i = Q~ one has d^ = 0, d^ ~ 0.68. 



B. Thermodynamic properties of hadronic phase 

Below we present numerical results for a nonstrange {ng = 0), isospin-symmetric hadronic 
matter. Figure [T] shows pressure P as a function of T for different values of //. One can 
see that the excluded volume corrections (EVC) reduce the pressure as compared to the 
calculation with v — 0. The difference between two calculations increases with raising T 
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TABLE I: Characteristics of hadronic species included in the calculation a . 
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a Values with asterisk correspond to antihypcrons. 
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FIG. 1: Pressure of nonstrange hadronic system as a function of temperature at different values 
of the baryon chemical potential \x . Thick and thin lines correspond to the excluded volume values 
v = 1 fm 3 and v = 0, respectively. 

and fi. At \i > 1 GeV it is large for any T. Pressure isotherms as functions of fi are shown 
in Fig. EL In accordance with Eq. (|T0l) . smaller derivatives of the isotherms imply lower 
baryon densities at v ^ as compared with the case v = . This is indeed seen in Fig. [31 
As expected, the saturating behavior n —* 1/v is predicted at fi — ► oo . 

To demonstrate sensitivity of the hadronic EOS to the strange chemical potential in 
Fig. H] we compare isotherms of pressure calculated assuming either n,s = or /ig = (in 
the latter case the strangeness neutrality is not guaranteed) . One can see that deviations 
caused by nonzero fis increase with /z, but they are rather small, especially at low T. 
Comparison of /i^-isotherms as functions of fi, calculated with and without EVC, is shown 
in Fig. [5l For v ^ the model predicts almost linear relation between fi$ and \i with the 
slope of about 1/2. The deviations from the ideal gas of point-like hadrons (v = 0) become 
important at fi > 1 GeV. 

In connection with nonzero fis one has to study the possibility of the if-meson con- 
densation in a dense hadronic medium. This problem is under investigation for more 
than two decades 23j. As well known, the Bose condensation of the i-th mesonic species 
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FIG. 2: Pressure isotherms for nonstrange hadronic system as functions of the baryon chemical 
potential fi . Thick and thin lines correspond to values v = 1 fm 3 and v = 0, respectively. 
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FIG. 3: Same as Fig. [2] but for isotherms of baryon density. 



becomes possible if the effective chemical potential (HI) exceeds the meson mass rrii 60]. 



Within the considered model one may expect positive chemical potentials /ij only for 
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FIG. 4: Pressure isotherms of hadronic system. Thick and thin lines are calculated assuming 
ns = and /U5 = 0, respectively. All curves correspond to v = 1 fm 3 . 
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FIG. 5: Isotherms of fis for nonstrange hadronic system: thin and thick lines correspond to v = 
and v = 1 fm 3 , respectively. 

strange mesons with Si = 1. In this respect the most "dangerous" are K + mesons with 
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FIG. 6: The effective chemical potential of kaons as a function of // at different temperatures T. 
Thin and thick lines are calculated assuming v = and v = 1 fm 3 , respectively. The dashed-dotted 
line shows the threshold of the if + -meson condensation. 

mass rriK+ — 494 MeV. As one can see from Fig. El in the case v — the Bose-condensation 
of K + mesons would be possible at /i > 1 GeV. However, in the calculation with high 
enough v, the effective chemical potential JIk does not exceed the condensation threshold 
at any /i and T. Thus, the EVC result in total suppression of the i\" + -condensation in 
chemically equilibrated hadronic matter with n$ = . It would be interesting to repeat this 
analysis for nonzero ns expected in stellar environments. 



C. Hadron multiplicity ratios in heavy ion collisions 



In this section we discuss an interesting behavior of the K + /n + and A/ V mu 
ratios observed in relativistic heavy-ion collisions by the NA49 Collaboration p], 



tiplicity 
SI. The 



experimental data show a peak ("horn") at low SPS energies. Our consideration follows 



closely the thermal model used in Refs. [2 



It is assumed that production of secon- 



dary hadrons in a nuclear collision can be described as the emission from a statistically 
equilibrated volume of hadronic matter characterized by certain temperature and baryon 
chemical potential. Within this model the kaon to pion multiplicity ratio equals n* K /n^ 
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TABLE II: Temperature and chemical potentials obtained by thermal fit of hadron ratios in central 
Au+Au and Pb+Pb collisions at different cm. bombarding energies y / ijv r /v. 



y/s NN , GeV 


T, MeV 


H, MeV 


»S, MeV 


2.70 


67 


735 


122.5 


3.32 


85 


668 


115.7 


3.84 


97 


621 


111.7 


4.30 


106 


584 


108.5 


4.85 


117 


545 


107.6 


6.41 


132 


460 


94.9 


7.74 


141 


405 


86.4 


8.87 


147 


368 


80.9 


12.4 


155 


287 


64.0 


17.3 


159 


219 


48.5 


62.4 


164 


69 


15.1 


130 


164 


34 


7.4 


200 


164 


22 


4.8 


5500 


164 


0.8 


0.2 



where n* K and n* are, respectively, the equilibrium densities of kaons and pions including 
those produced in decays of hadronic resonances. We apply the following relations 



n 



K 



m. 



(17) 



Here rij = n.j (/i, /is, T) is the partial density of i-th hadrons at fixed /i, fi$ an d T calcu- 
lated by Eq. (fl3l) . As before, /is is determined from the condition of strangeness neutral- 
ity (jHJ). The numerical values of dj and df are listed in Table [B By definition, these factors 
equal unity for directly produced pions and kaons, and they are zero for stable hadrons 
with % ^ tt,K. In our calculations we use the chemical freeze-out parameters //, T deter- 
mined from thermal fits of hadronic ratios measured at the AGS, SPS and RHIC energies. 
Namely, we apply the parametrizations of /i, T as functions of the bombarding energy ^/snn 
given in Ref . 29( . In this way we get the values listed in Table [TTJ The last column of the 
table gives fj,$ = HsifaT) calculated by using Eqs. (JHJ) — (JSD , (IS])- (fTT| for the case v = lfm 3 . 
The last line is obtained by extrapolating the parametrizations of /i, T to the LHC energy. 

In Fig. [7] we show the densities n* , n* K as well as the baryon density (at freeze-out stage), 
calculated by using the T, /i, /is values from Table ILT1 According to this calculation, the max- 
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FIG. 7: Total densities of pions (thick solid line) and kaons (thick dashed line) at chemical freeze- 
out in central heavy-ion collisions as functions of cm. bombarding energy (v = lfm 3 ). Thin lines 
give the contributions to pion and kaon densities from resonance decays. The dashed-dotted line 
shows the net baryon density. 
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FIG. 8: The K + /tt + multiplicity ratio in central nuclear collisions as a 
energy (v = lfm 3 ). Experimental data are taken from Refs. [7I. I^. I30 . 



function of cm. bombarding 
31]- 



imal baryon density of about 0.5 n is reached at ^snn — 9 GeV (i?i a b = 40 AGeV). One can 
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see that the pion density saturates at ^/snn > 40 GeV on the level of about 0.3 fm~ 3 . On the 
other hand, the maximal kaon density ~ 0.04 fm -3 is reached already at ^snn ~ 10 GeV. 
It is interesting to note that at high energies the contributions of resonance decays reach 
about 70% and 50% in the case of pions and kaons, respectively. 

In Fig. [H]we compare our results for the K + /ir + ratio 6l|] with experimental data obtained 
at the AGS [30], SPS [7, 8] and RHIC jsi| energies. The solid line represents our standard 
calculation, where the broad /o(600) resonance (a meson) is disregarded. One can see that 
our model is able to reproduce qualitatively the observed Kj 7r ratios but the agreement with 
data is not perfect. In particular, the experimental data exhibit a much sharper peak at 
y/SNN — 7.74 GeV. Choosing different values of v gives only slight modifications and does 
not improve the shape of the Kj 7r excitation function. This is easy to understand since the 
EVC are practically cancelled out in hadron multiplicity ratios. 
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FIG. 9: The A/ir~ multiplicity ratio in central nuclear collisions as a function of cm. bombarding 
energy (v = 1 fm 3 ). The dashed lines are calculated with inclusion of additional mean- field potential 
of A hyperons. Experimental data are taken from Refs. [?], ^, I30I 
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Recently, the authors of Ref. [29(] have noticed that these data can be better reproduced 
by the thermal model which takes into account the contribution of the /o(600) resonance. To 
verify this observation, we have made additional calculations including this meson. These 
calculations are done in the zero-width approximation assuming different masses of a me- 
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son, and using the values g a = 1, d£ = 2, = 0. As seen in Fig. [8] visible changes in 
the K + /n + ratio take place at ^snn > 5 GeV. However, even at smallest m a ~ 0.4 GeV, 
due to increased pion multiplicity, this ratio drops only by about 10%. At the same time, the 
observed values in the peak region are underestimated stronger in this case. As expected, 
the shifts of the K + /n + ratio become smaller for larger m a . On the basis of this analysis, 
we conclude that the inclusion of the a meson does not solve the problem. According to 



Ref. (34j, to obtain a better fit of data one should introduce an additional parameter 75 < 1 
which is responsible for the suppression of strange hadron yields compared to the equilibrium 
model predictions. 

We have also calculated the excitation function of the A/tt~ ratio. The results are shown 
in Fig. One can see that the model can qualitatively reproduce the experimental data, al- 
though the peak height is somewhat underestimated. Our analysis shows that the agreement 
with the data can be improved if one takes into account an additional mean-field potential 



of A hyperons, U\ = ^U(n) 62|, where U(n) is the corresponding Skyrme-like potential for 
nucleons introduced in Sect. IV A[ 

Finally, we would like to emphasize that a nonmonotonous behavior of the K/tt and A/ V 
excitation functions is obtained without any reference to the QGP. Within our model this 
anomaly is mainly related to the strangeness neutrality condition which in turn requires 
nonzero (is in the hadronic phase. In the Boltzmann approximation the K/tt and A/ it 
ratios are approximately proportional to i/k and y\ defined as 

y K = exp [((i s - m K + m % )/T] , y\ = exp [(/i - (i s - m A + m n )/T] . (18) 

According to our calculations these functions have maxima at ^/snn — 9 and 5 GeV, 
respectively. The appearance of these maxima follows from a nonmonotonous behavior of 
the strangeness fugacity exp((i s /T) along the chemical freeze-out line [6^]. As demonstrated 
above, this behavior agrees qualitatively with experimental data. The inclusion of resonance 
decays makes the K/tt and A/n peaks more pronounced. 

III. QUARK GLUON PHASE WITHIN THE BAG MODEL 

To calculate thermodynamic properties of a baryon-rich QGP we use a simple bag model 
with perturbative corrections of the order a s . Gluons (i = g) and light quarks (i = q, q) 
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are considered as massless point-like particles, but for strange quarks (i = s,s) we intro- 
duce a nonzero mass m s . Nonperturbative effects are introduced via the bag constant B . 
In the chemically equilibrated QGP, one can again express the chemical potential of the 
2-th particle by Eq. (jSj), where fi and fis are now the baryon and strange chemical poten- 
tials of the quark-gluon system. As a result, we get the relations \i g = 0, \i q = —fig = /i/3, 
fi s = -lis = A*/3 - up. 



Following Ref. 



351 ] we take into account perturbative corrections by introducing additional 



constant factors 1— £ and 1 — 0. 8£ into the kinetic pressure of quarks and gluons, respectively. 
Here £ ~ a s is the model parameter which we fix by the comparison with lattice data. Within 
such an approach one gets the following expression for pressure of the QGP: 



oo 



1 9 2 V 90 ' V 18 3247r 2 y 7T 2 y v ; le^ + 1 



e t + 1 J 

Here iV 9 = 16(1 — 0.8£) is the effective number of gluons and Nf = 2(1 — £) is the 
effective number of light flavors. The third term in Eq. (fTDl gives the contribution of 
strange quarks and antiquarks. Corresponding expressions for densities of baryon charge n, 
strangeness ns, energy e and entropy s can be obtained from Eq. (fl9|) by using the ther- 
modynamic relations ([7]), ffT4l) . One can easily see that the strangeness density of the QGP, 
n s — n s — n s = 9P/ dfis > vanishes if [i s = us ~ A*/3 = . 

Below we fix m s = 150 MeV and take the same value B = 344MeV/fm 3 as in our previ- 



ous fluid-dynamical calculations (EOS-I) 36|, |37| of heavy-ion collisions at RHIC energies. 



We have calculated the thermodynamic functions of the QGP for different values of the 
parameter £. Motivated by the comparison with lattice data (see Fig. [TTj) . we further use 
the value £ = 0.2 in our calculations. The advantage of the bag model is that it can be used 
in the region of nonzero chemical potential, which is still not accessible by lattice calcula- 
tions. In Fig. [10] we show pressure isotherms of strange (us — 0) and nonstrange (fis = /x/3) 
QGP. Similarly to the hadronic matter (see Fig. Hj), in the case of nonzero fis t ne pressure 
is reduced compared to the calculation with fi S = 0. 
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FIG. 10: Same as Fig. HI but for pressure isotherms of the quark-gluon plasma. Notice the 
appearance of negative pressure values at low T and [i . 

IV. THE DECONFINEMENT PHASE TRANSITION AND PROPERTIES OF 
MIXED PHASE 

A. Phase equilibrium conditions 

Below we use our phenomenological approach, where the hadronic and quark-gluon phases 
are described by two different models, to study the possibility of a first order deconfinement 
phase transition. If such a transition exists, then only three types of equilibrium states are 
possible at a given point in the (//, f-is,T) space: the hadronic phase (HP), the quark-gluon 



phase (QP) and the mixed p 
which has a higher pressure 



lase (MP). Outside the MP region the stable phase is the one 



121 ]. The Gibbs condition of equilibrium between the domains 



of different phases in the MP can be written as 



Ph (fJ>, fis, T) = Pq (/!, us, T) 



(20) 



Here and below we mark thermodynamic functions of the hadronic and quark-gluon phases 
by indices H and Q, respectively. The expressions for pressure Ph and Pq are given in 
Sect. EDS and EH 

In a baryon-rich mixed phase the condition of zero total strangeness leads to the 
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strangeness-antistrangeness separation phenomenon 20]. At nonzero fi and T, the equa- 
tions ush — n SQ — may hold simultaneously with Eq. fl20|) only in a single point of 
the (fi, [is, T) space. But in general, the strangeness numbers of coexisting domains are 
nonzero and compensate each other only on average over the total volume of the system. 
Let A denotes the volume fraction of hadronic domains in the MP: 

^v^M- < 21 > 

where Vh and Vq are the total spatial volumes occupied by hadrons and the QGP, respec- 
tively. Then the condition of zero total strangeness may be written as 

n s = XnsH + (l-X)n S Q = 0. (22) 

This condition holds if n SH and usq have different signs. At fixed [i, A the phase transition 
temperature T c is determined by solving Eqs. ( |20l) . ( |22l) . The phase transition region in 



the \i — T plane is not a line, but a strip 21[ T = T c (/i, A) where T c is decreasing function 



of A . The lines T = T c (/i, 1) and T = T c (/z, 0) give the hadronic and quark-gluon boundaries 
of the MP. Our calculations show (see Fig. [15]) that the strip's width decreases with increasing 
excluded volume v. 

The energy, baryon and entropy densities in the MP are calculated by using the first 
equality of (1221) with the replacement of ns by e, n and s , respectively. At given /i, T we use 
the values of fis, A obtained by solving Eqs. (12"U|) . (12"21 . 



B. Results for baryon-free matter 

In this section we present our results for the baryon-free system (n — jj,s — 0), which 
can be compared with lattice calculations. Table II I II gives critical temperatures T c obtained 
for different values of £. In Fig. [TT] we compare the model predictions for e/T 4 with the 

TABLE III: Critical temperature of the deconfinement phase transition in the baryon-free matter 
(v = 1 fm 3 , B = 344MeV/fm 3 ) . 
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lattice calculations from Ref. [lj. One can clearly see a pronounced peak just above T c 
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FIG. 11: Scaled energy density of baryon-free matter as a function of scaled temperature at 
different values of £. All results are obtained with v = lfm 3 and B = 344MeV/fm 3 . T c is the 
critical temperature of the deconfinement phase transition (see Table HIB . Open and closed circles 
show the lattice data [1] with the number of quark flavors Nf = 2 and 3 , respectively. 



which is not present in the lattice data. This is a well-known artefact of the bag model (see 
e.g. {37]) ■ At high temperatures a good agreement with the lattice data can be achieved by 
choosing £ ~ 0.2 . Based on these results, we use £ = 0.2 in the following calculations. 

A more detailed information about the EOS of baryon-free systems is given in Figs. H2HHL 
Figure [12] represents the graphic solutions of Eq. (120!) for zero and nonzero v. In the 
case v = 1 fm 3 the curves Ph{T) and Pq(T) intersect in a single critical point C . On the 
other hand, in the calculation with v — two crossing points C and C" are present. The 
second point corresponds to the transition QGP — > HP at the temperature T = Tc» > Tc- 
Such unusual behavior, not supported by lattice data, is explained by a too steep rise of 
pressure of point-like hadrons with increasing number of species at high temperatures. 

To perform fluid-dynamical modelling of heavy-ion collisions, one should know pressure 
as a function of the energy density e and net-baryon density n . Such a representation of the 
EOS of the baryon-free matter is shown in Fig. [131 In this figure the deconfinement phase 
transition shows up as a horizontal line connecting the HP and QP. This line corresponds to 
the MP states. The parameters of the phase transition, calculated for v — 1 fm 3 , are close to 
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FIG. 12: Pressure of baryon-free matter as a function of temperature. The solid and dotted lines 
show pressure of the HP calculated with v = 1 fm 3 and v = 0, respectively. The dashed line shows 
pressure of the QGP calculated with B = 344MeV/fm 3 and £ = 0.2. 



those in the equation of state EOS-I used in Refs. |36|, |37|. The dashed line shows metastable 
states of the HP extended into the region of large e. At such energy densities, a much harder 
EOS of the HP is predicted compared to the calculation with v = (shown by the dotted 
line). On the other hand, both calculations giv e similar results at low e. In fluid-dynamical 



simulations of heavy-ion collisions 



37 



381 ] two scenarios, with and without the phase 
transitions, are often compared to check the sensitivity of observables to the EOS. However, 
from comparing the dashed and dotted lines in Fig. [13j it is evident that such analysis could 
be misleading if the role of repulsion in a dense hadronic system is ignored. 

The sound velocity c s is an important characteristic of the EOS which gives the speed of 
small perturbations of matter in its local rest frame. Within the ideal hydrodynamics the 
sound velocity squared is equal to |l_2j 



dP\ fdP\ n (dP\ 
fe)a \ de )n w\dn) t 



(23) 



where a = s/n is the entropy per baryon and w = e + P is the enthalpy density. In the 
second equality we have applied the thermodynamic relation (n s = 0) : 



da 



1 

nT 



' de 



w 



n 



dn) . 



(24) 
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FIG. 13: Pressure of baryon-free matter as a function of energy density. Dashed and dotted lines 
show pressure of a hadronic system calculated with v = 1 fm 3 and v = 0, respectively. The solid 
line shows pressure calculated with inclusion of EVC and phase transition effects. Above the dot 
the HP has superluminal sound velocities. 
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FIG. 14: Same as Fig. [131 but for sound velocity squared. 
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In the case of baryon-free matter (n = 0) is equal to the slope of pressure as a function 
of e. 

The causality condition c s < 1 should be fulfilled in any model consistent with relativistic 
kinematics. However, a simple Van der Waals approach does not guarantee this property, as 



has been already mentioned in Refs. [16|, ll7|]. In a gas of hard spheres the acausal behavior 
is associated with an implicit assumption that the signal propagates instantaneously over 
the sphere extension. As will be shown below, the causality condition c s < 1 is indeed 
violated in the HP, but only at rather high baryon densities and small temperatures. For 
illustration, Fig. [UP shows c 2 s for baryon-free matter. One can see that c 2 s is close to 1/3 
in the QP, vanishes in the MP and equals 0.2 — 0.25 (at v = 1 fm 3 ) in the HP. The values 
c s > 1 are reached only in metastable hadronic states with e > 1 GeV/fm 3 which are not 
realized in the equilibrated matter. 

C. Phase diagram in \i — T plane 
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FIG. 15: Phase diagrams in the fj, — T plane calculated for different v. Thick and thin lines 
correspond to hadronic and quark-gluon boundaries of the MP, respectively. 

Using the Gibbs condition fT20l) we have calculated the phase diagram of matter in 
the fi — T plane. The results for different values of v are given in Fig. [15j For each v 
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we show two boundaries of the MP corresponding to the conditions nsu = (A = 1) 
or tisq = (A = 0). One can see that the width of the MP region decreases with in- 
creasing v. The two boundaries are practically undistinguishable at v > 0.5 fm 3 . With 
decreasing v the HP occupies larger and larger domain of the fi — T plane. At v < 0.1 fm 3 
the phase boundary exhibits a back bending at small chemical potentials, < fx* ~ 0.8 GeV. 
This means that at fixed /i < /i* and increasing T three phase transitions appear: the first 
and third ones from the HP to the QGP, and the second intermediate transition, from the 
QGP to the HP. The third point goes to infinity at v — > 0. This means that the HP is ther- 
modynamically more stable than the QP at asymptotically high temperatures. At v = 
only a small domain of the QGP (its boundary is shown by the dotted line) remains in 
the jJL — T plane. This is certainly an unphysical behavior, which clearly shows inconsistency 
of the hadron resonance gas model with point-like hadrons. 
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FIG. 16: Strangeness densities of the hadronic (thick lines) and quark-gluon (thin lines) compo- 
nents in the MP region of the phase diagram as functions of the baryon chemical potential. The 
results are shown for several values of the volume fraction A . For A = (no hadrons) hsq = . 



Figure [16] shows the strangeness number densities ush and usq across the MP (v = 1 fm 3 ) . 
One can see that indeed, n SH and n S Q have different signs (see Sect. IIV Ap . but their 
absolute values do not exceed 0.1 fm -3 . At \i > 0.8 GeV (anti) strangeness in hadronic 
domains is carried mostly by hyperons. According to Fig. [T71 at large chemical potentials 
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FIG. 17: Strangeness density of the quark-gluon domains in the MP (A = 0.5) at different values 
of the excluded volume v . 

the strangeness density usq is rather sensitive to v . 



D. Adiabatic trajectories 

Further on we consider in more details the results for v — 1 fm 3 . In Figs. [T8H221 we 
show the adiabatic trajectories (" adiabates" ]i.e. sets of states with equal entropy per 
baryon, S / B = a = const. As well-known [39y, in the ideal hydrodynamics the entropy 
per baryon is conserved in a given fluid element. Therefore, these trajectories contain an 
important information about the conditions which can be realized in heavy-ion collisions. 
Figure [T8l shows the phase diagram and corresponding adiabates in the [i — T plane. In such 
a representation the adiabates with S/B = and S/B = oo are given, respectively, by the 
horizontal (T = 0) and vertical (fi = 0) axes. Note that at T — ► all adiabates with finite /i 
end in the point fi ~ mjv, where = 939 MeV is the nucleon mass. Moreover, these 
adiabates have a zigzag-type behavior characteristic for a first order phase transition. This 
means that along the adiabatic trajectory the temperature grows when the system enters the 
coexistence region. In other words, the temperature of the HP at A = 1 is higher than the 
temperature of the QP at A = 0. Such a picture differs from the predictions of the linear a 
model and the Nambu-Jona-Lasinio model 0, |4l| for the chiral first order phase transition. 
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FIG. 18: Lines of constant entropy per baryon in the \i — T plane [v = 1 fm 3 ) . The solid curve 
represents the phase transition line for A = 0. Full dots corresponds to the /i, T values obtained 
from thermal fits of hadron yields 4CJ] observed in central Au+Au and Pb+Pb collisions at dif- 
ferent bombarding energies. The region between the dotted and thin solid lines contains states 
with c s > 1 . 



There temperature drops in the MP. This difference may be related to the fact that the 
present model includes massless gluons which are completely ignored in Refs. 0,0. 



Full dots in Fig. [18] show states of the chemical freeze-out [42J in central collisions 
of heavy nuclei with cm. bombarding energies from GSI (y / J/v r /v — 2.3 GeV) to RHIC 
{\/snn — 200 GeV). "Experimental" values of T and /i have been found 40J from thermal 
fits of hadron multiplicity ratios observed in such collisions. 

A word of caution is in place here. As has been already mentioned, our model contains 
superluminal sound velocities in the HP. These states are situated in the region between 
the dotted and thin solid lines in Fig. [181 One can see that such states correspond to 
values S/B < 5. On the other hand, it is known from hydrodynamical simulations (see 
e.g. 43J) that typical conditions realized in heavy-ion collisions at energies > 10 GeV 
correspond to entropy per baryon S/B > 10. Therefore, an adiabatically expanding system 
of particles produced in ultrarelativistic nuclear collisions does not enter the region of c s > 1. 
However, such states can be reached in compact stars. Based on these results, we conclude 
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FIG. 19: Boundaries of different phases and adiabatic trajectories in the n — T plane (v = 1 fm 3 ) . 
The shaded area shows the mixed phase region. The hadronic states on the right from the dotted 
line have sound velocities c s > 1. 
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FIG. 20: Same as Fig. [19] but in the n — e plane, 
that the present model should be considerably modified to study properties of the HP and 
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the deconfinement phase transition at high baryon densities (see the discussion in Ref. [9(). 
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FIG. 21: Adiabatic trajectories in the e — P plane (y = 1 fm 3 ) . 



3.0 



o 



1.0 - 



0.5 - 



V 


I 


1 I 1 I 


1 






v=1 fm 3 




II 




— S/B = 
5 

30 

oc 




■ // 
- ' k 

i! 














1*~-T- 


, — u , ! . 





e(GeV/fm d ) 

FIG. 22: Sound velocities squared in adiabatic processes with different S/B as functions of energy 
density (v = 1 fm 3 ) . 



In Figs. [T914201 we show boundaries of different phases in the n — T and n — e planes. 
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Again, one can see that adiabates with S/B > 10 do not enter the "dangerous" region 
with superluminal sound velocities. The solid line in Fig. [20] corresponds to the limiting 
case of T = 0. States below this line can not exist. Figure |2"T1 shows pressure adiabates 
versus the energy density. According to first equality of (j2"5|) . c 2 a - values are given by slopes 
of adiabates in this representation. Note, that pressure in the MP is constant along the 
lines fi = const which differ from the lines cr = const in the n — e plane. The direct 
calculation of function of e gives the results shown in Fig. [22j One can see that c 2 s 

values especially large at zero temperature, for states near the MP boundary. Generally, 
the sound velocities are nonzero in the MP, however, they are much smaller than maximal 
values of c s in the HP. The nontrivial two-peak structure of the adiabate with S/B = 
appears due to the contribution of baryonic resonances (A, N* . . .). However, this behavior 
takes place only in the unphysical region with c s > 1. 



V. INTRODUCTION OF THE MEAN FIELD 

A. Skyrme-like parametrization of the mean-field potential 

The EOS obtained in preceding sections takes into account only short-range repulsive 
interactions of hadrons and completely ignores the intermediate-range attractive interaction 
between baryons. Due to this reason, this EOS does not describe the saturation property of 
isospin-symmetric nuclear matter at low temperatures. It is well-known that such attractive 
interaction leads to the liquid-gas (LG) phase transition in nuclear matter at tempera- 
tures T < lOMeV and chemical potentials /x ~ /i = — E B , where E B ~ 16 MeV is 
the binding energy of cold nuclear matter. This phase transition manifests itself as the 



multifragmentation phenomenon in intermediate-energy nuclear reactions [441 ]. 

To account for the mean-field effects, we introduce an effective potential U = U (n) which 
depends only on the baryon density n and does not depend on momenta of interacting 
baryons 64|. Then the baryon's single-particle energy can be obtained simply by adding U(ri) 
to the kinetic energy. In this case the partition function of the hadronic system can be 
calculated analytically ]16]. Again, finite sizes of hadrons are taken into account by the 
volume reduction (TfJ). As the result, the following formulae for thermodynamic functions of 
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the HP can be written 



A« = m + U(n) , 
P = P x + P / (n) 
e = e K + e f (ri) . 



(25) 
(26) 
(27) 



The "field" contributions (marked by index /) to the energy density and pressure are found as 



€f(n) = nU(n) — Pf(n) = / dn\U{ni) 



(28) 



The kinetic terms (marked by index K) in Eqs. (I2"5"l) - (l2"7j) are functions of hk,Hs,T cal- 
culated by using formulae of Sect. Ill Al with the replacement /i — > fi K , P — > P K , e — > e K . 
The resulting expressions for P K and ex are obtained from Eqs. ((3j)-((SD and (iTST) . respec- 
tively. One can show that the total pressure P = P((jl, Hs,T) satisfies the relation ([7|). The 
densities s,n and ris are given by Eqs. (l9l- (fTTl) where now 



Hi = Bifx K + Siiis - vP< 



K 



(29) 



In a spirit of the Skyrme approach [45] we parametrize the mean-field potential in the 
form 

u{n) = - a (£} + e(£j , ( 30) 

where n is the saturation density of nuclear matter and a, j3, 7 are density-independent 
parameters. In the following we fix 7 to a commonly used value 7/ 6 and choose the remaining 
parameters from the requirements P = 0, e/n = /io = 923 MeV at n = no, T = 0. The 
values of a, (3 as well as the incompressibility modulus K = 9dP/dn at the saturation 
point, calculated for different choices of v, are given in Table HVl 

TABLE IV: Parameters of the mean-field potential and the incompressibility modulus of equilib- 
rium nuclear matter for different values of the excluded volume. 



v, fm 3 





1 


2 


a, MeV 


352 


334 


297 


0, MeV 


301 


277 


230 


K, MeV 


200 


214 


264 
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FIG. 23: Pressure of the baryonic system at T = as a function of /i for different values of the 
excluded volume v. Thick and thin lines correspond to the cases U = U(n) and U = 0, respectively. 
The part DA of the thick solid curve represents unstable states for v = 1 fm 3 . 




n (fm" 3 ) 

FIG. 24: Same as Fig. 1231 but for pressure as a function of the baryon density n. 

Figures I23ti25l illustrate the properties of cold nuclear matter as predicted by this model. 
Figure [23] shows pressure as a function of /i for several values of the excluded volume v. 
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FIG. 25: Same as Fig. [2H but for energy per baryon. 

Point B marks the saturation point of cold nuclear matter. The mean-field is responsible 
for the appearance of several branches of pressure at fi < fiA — m N • In the case v — 1 fm 3 
the branches BA and DB describe the metastable states, while the branch DA corresponds 
to the unstable states (see below). This behavior differs qualitatively from the case U = 
where only one branch is present at fi ^ [ia- Figures l2~4TT2"5l show pressure and energy per 
baryon as functions of n. By construction, the saturation point has the same position for 
all considered values of v . 

B. The liquid-gas phase transition 

Now let us consider properties of hadronic matter at nonzero temperatures. Figure [261 
shows isotherms in the \i — P plane calculated with v = lfm 3 . In this case the existence 
of the critical point is predicted at the temperature T c ~ 14.925 MeV and the chemical 
potential /i c ~ 909 MeV (see Table IVl) . The isotherms with T <T C contain three branches 
of pressure, similarly to the case T = discussed above. The phase transition points 
correspond to intersections of metastable parts of the isotherms. The branches with negative 
curvature correspond to unstable (spinodal) states. Small density perturbations will grow 
exponentially in this region. This is a well-known spinodal instability which leads to the 
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FIG. 26: Isotherms in the n — P plane (v = 1 fm 3 ) . 
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FIG. 27: Pressure isotherms as functions of baryon density (u = 1 fm 3 ) . The shaded area shows 
the mixed phase region of the LG phase transition. The dotted line with maximum at point C is 
the boundary of the spinodal instability domain. 



separation of matter into dense and dilute domains 46(, characteristic of the LG phase 
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transition. We would like to note here that the LG mixed phase can be represented by the 
ensemble of nuclear fragments of different sizes. This ensemble can be well described by 
the statistical multifragmentation model which also includes the EVC. The anomaly 
in the caloric curve, associated with the LG phase transition, has been indeed observed in 
intermediate-energy heavy-ion collisions by the ALADIN collaboration 48]. 

Figure [271 represents the isotherms in the n — P plane. The shaded region corresponds to 
the mixed phase, where the hadronic gas (nucleons and light clusters) coexists with nuclear 
fragments (droplets of liquid). Point C is the critical point of the LG phase transition. The 
parameters of the critical point for different values of v are given in Table |Vj One can see 



TABLE V: Characteristics of the critical point of the LG phase transition at different values of the 
excluded volume. 



v, fm 3 





1 


2 


T c , MeV 


15.37 


14.93 


12.67 


Mc MeV 


907 


909 


915 


n c , fm~ 3 


0.048 


0.048 


0.039 


P c , MeV/fm 3 


0.198 


0.195 


0.139 



that their sensitivity to the excluded volume v is rather weak. This is also seen in Fig. 1281 
where we compare spinodals and the phase transition lines for v = and v = lfm 3 . On 
the basis of these results we conclude that characteristics of the LG phase transition in our 
model are similar to predictions of other authors (see e.g. 491]). 



C. Phase diagram with two phase transitions 

Finally we present results for the full model where both the LG and deconfinement phase 
transitions are included. In the case T = 0, applying Eqs. (|25]) - (|50]) for a broad interval 
of baryon densities we get the results shown in Fig. [291 According to our calculations, 
in the region of the deconfinement phase transition the corrections due to the mean-field 
interaction are rather small, of the order of 5% . Note that the mean-field potential leads to 
vanishing pressure for zero temperature states with n < n (see Fig. |2"7j) . 

Figures EDHSIl represent the full phase diagram in the fi—T and n—T planes for v = 1 fm 3 . 
Compared to the calculation with U = 0, a new phase transition line starting at /i = /io 
appears in Fig. [301 As seen in Fig. SB the LG mixed phase (shown by shading) occupies a 
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0.95 



H (GeV) 

FIG. 28: Phase diagram of the LG phase transition in the [i — T plane (the solid line). The 
dashed lines show boundaries of the spinodal region. Thin and thick lines correspond to v = and 
v = lfm 3 , respectively. 

relatively small region of the n — T plane. According to Fig. [31] the borders of the quark- 
hadronic mixed phase are shifted only slightly due to the mean-field effects. The same 
conclusion is valid for the mixed phase boundaries and the T = line in the n — e plane 
(cf. Fig. [20]). 

VI. SUMMARY AND DISCUSSION 

In this paper we have investigated the EOS of strongly interacting matter in a phe- 
nomenological model taking into account the excluded volume effects in the hadronic phase. 
The quark-gluon phase is described by the bag model with the lowest order perturbative 
corrections in kinetic terms. Within this approach one can get only a first order phase 
transition. The sensitivity of the phase diagram to the excluded volume v has been inves- 
tigated in details. At small values of v the phase diagram has unphysical behavior since 
the hadronic phase becomes preferable at high temperatures. Only at v > 1 fm 3 the criti- 
cal temperature of the deconflnement phase transition monotonously decreases with baryon 
chemical potential /x. Generally, the strength of the phase transition increases with increas- 
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FIG. 29: Pressure as a function of baryon density at T = 0. The dashed and solid lines are 
calculated with and without the mean-field effects, respectively (v = lfm 3 ). Note change of the 
vertical scale compared to Figs. [Ml [27J 




FIG. 30: Phase diagram of strongly interacting matter in the [i — T plane (v = 1 fm 3 ) . The dashed 
and solid lines are obtained from calculations with and without the mean-field effects, respectively. 

ing /i. Typical values of the baryon density in the mixed phase are in the range (5 — 10) no 
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FIG. 31: Same as Fig. [30] but in the n — T plane. Shading shows the region of the LG phase 
transition. 

at T ~ 100 MeV. Such a strong phase transition should certainly lead to clear observable 
signals, e.g. the formation of quark-gluon droplets at final stages of the relativistic nuclear 



collisions 
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5l|- 



The condition of zero net strangeness has been imposed in our calculations. It leads 
to almost linear relation between the strange and baryon chemical potentials. It is shown 
that the Bose condensation of K + mesons is not possible at realistic values of v. We have 
demonstrated the possibility of the strangeness-antistrangeness separation in t he mixed 
phase. This effect may help to produce clusters of strange matter like strangelets 
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MEMOS 



531 ] or 



54J ] in relativistic heavy-ion collisions. We have calculated the adiabatic trajecto- 



ries (S/B = const) and found that they have a zigzag-like shape in the mixed phase region. 
According to our analysis, sometimes the model predicts acausal states (c s > 1), but they 
lie outside the region reachable in ultrarelativistic heavy-ion collisions. By implementing 
the mean-field potential for baryons we have described simultaneously the liquid-gas and 
deconfmement phase transitions. However, properties of the quark-hadron mixed phase are 
only slightly influenced by the mean-field effects. 

In the future we are going to generalize this model by introducing different excluded 
volumes for different hadronic species. Some attempts in this direction have been already 
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made in Refs 



55] . Our present approach does not take into account modifications of 



particle properties in a dense medium, although hadronic masses and radii may significantly 
change as compared to their vacuum values. This problem should be also studied in the 
future. 
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APPENDIX A: THERMODYNAMIC FUNCTIONS OF HADRONIC SYSTEM 



In the case of fermions (Bi = ±1) one can represent the integrals in Eq. (fl6|) in the 
form [12| which makes easier their numerical calculation at low temperatures: 
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Here pp = a/ fif — mf , (x) = \ (1 + sgnx) and ijj (x) is defined as 
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The integrals in the first term of Eq. (1A.1I) vanish at T — > . They were calculated numeri- 
cally using the Newton-Cotes method. 

In the case of mesons (Bi = 0) we calculate integrals in Eq. ffl6l) by representing them as 
series of modified Bessel functions: 



Pi 



9i™>i \^ I Hi 1 

"2^ exp It" 



i=i 



mjK 2 
x 2 

X 



(A.3) 



37 



where x = rriil/T and K n = K n (x) is the MacDonald function of the n— th order. These 
series converge if /I, < rrii (see first footnote on page [6]). 
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As claimed in Ref . [3] , crystallization into a solid state should occur at densities near the close 

packing limit n cp = (4\/2r^) _1 ~ 3v~ l . Our estimates show, however, that a dense hadronic 

system should earlier undergo the transition into a quark-gluon phase. 

In Eq. (|16p we assume that the condition of Bose-condensation /Ij > nii is not satisfied (see 
below). 

Except of /o(600), a very similar set of hadrons has been used in the THERMUS thermal 
model 0. 

For the observed decay channels with unknown probabilities we assume equal branching ratios. 
Here we disregard possible in-medium modifications of meson properties. This question have 
been studied on the mean field level in Refs. 
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[61] Using the isospin symmetry we write K + /ir + = 1.5nt/rt*. 

[62] The factor 2/3 comes from a naive consideration based on the SU(3) flavor symmetry. 
[63] As was mentioned on page[U [i ~ \ /i^, so that \x — [is ~ /As- 

[64] For simplicity it is assumed that strange and nonstrange baryons have the same mean-field 
potentials. On the other hand, in the considered case of chemically equilibrated matter 
with zero strangeness, the medium-range interactions are not so important at high tem- 
peratures T > 100 MeV, when the abundances of hyperons, antibaryons and mesons be- 
come significant. 
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